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ABSTRACT 

We  show  that  the  effect  of  the  non-uniform  environment  on  the  static  screening  of 
the  Coulomb  potential  causes  a  major  change  in  the  calculated  surface  energy  of  simple 


metals. 


Accesjon  For 
NTIS  CRA&I 


Density  functional  theory*  (DFT)  needs,  for  practical  purposes,  approximations  to 
the  exchange-correlation  energy  functional  ExdPi  (**)].  Most  common  is  the  local  density 
approximation  (LDA)  which  is  valid  if  the  density  Pi(r)  varies  sufficiently  slowly  as  a 
function  of  r.  It  is  not  a  priori  clear  whether  this  condition  is  satisfied  in  a  metal  surface. 

As  an  alternative,  we  have  recently  (Ref.  2,  hereafter  called  paper  I)  developed  a 
variational  theory  for  inhomogeneous  Fermi  systems.  The  theory  has  been  shown  to  be  ac¬ 
curate  enough  to  provide,  in  the  homogeneous  limit,  resonable  agreement  with  the  known^ 
properties  of  the  bulk  system.  But  it  can  be  evaluated  without  further  assumptions  on 
the  slowness  of  the  density  variation.  We  have  applied  the  theory  to  the  calculation  of 
electron  densities,  surface  energies,  and  work  functions  of  metal  surfaces.  We  obtained 
surface  energies  of  idealized  jellium  systems  significantly  higher  than  the  ones  predicted  by 
density  functional  theory  (DFT)  in  the  local  density  approximation  (LDA)*.  This  letter 
reports  on  extensions  of  the  calculations  of  I,  and  clarifies  the  origin  of  the  larger  surface 
energies. 

The  variational  method  for  an  inhomogeneous  Fermi  system  starts  with  an  explicit 
ansatz  for  the  ground-state  wave  function 

i'l'o;  =  exp{i  ^u,(r,)  +  J  ^U2(r,,r;)}i<l>o)  (1) 

I  "  i<; 

<^n/  is  the  Slater  determinant  of  a  set  of  single-particle  orbitals  <pi[j)  =  0,(r;  )x(t),  (t-J  ~ 
1 . 4);  the  \(i)  are  the  spin-eigenfunctions. 


The  function  «i(r)  can  be  shown**  to  cancel  a  set  of  point-reducible  diagrams  if  an 
optimized  single-particle  basis  is  used,  and  U2(r, .r^)  is  determined  by  minimization  of  the 


FHNC  approximation  for  the  ground-state  energy, 

6E  6  {^o\H\^o) 

Su2  Su2  (^ol^o) 

The  single-particle  wave  functions  can  be  generated  by  a  one-body  equation  of  the 

form 

[-^V*  +  (/(r)|«.(r)  =  €,«(r)  (3) 

where  the  effective  one-body  potential  U{t)  is  determined  by  functional  minimization  of 
the  ground-state  energy^. 

=  0-  (4) 

This  procedure  is  equivalent  to  the  minimization  of  the  FHNC  ground-state  energy  with 
respect  to  the  one-body  density  in  the  space  of  all  densities  that  can  be  generated  by  a 
local  single-particle  potential.  In  that  sense  we  have  carried  out  the  DFT  program  without 
the  assumption  of  a  slowly  varying  density. 

The  theoretical  basis  of  I  is  the  Fermi-hypernetled  cliain  (FHNC/'/n)  theory  for  inho¬ 
mogeneous  systems'^.  Using  the  FHNC /  0  approximation,  we  have  derived  in  I  expressions 
for  the  correlation  energy  and  equations  determining  the  two-body  correlations  and  the 
optimal  single-particle  potential.  To  review  the  main  ingredients  of  the  theory,  we  need 
a  few  definitions:  The  convolution  products  of  two  functions  /l(ri,r2)  and  fl(ri,r2)  is 
written  as 

A  »  fli(ri,r2)  s  j  d^r3/l{ri,r3)i?(r3,r2).  (5) 

.''Xri .  r2)  is  the  static  form  factor  of  the  model  system  described  by  !^o))  inverses  like 


6E\U\ 

6U{r) 


2 


Sp  *(ri,r2)  are  understood  in  the  sense  of  the  convolution  product  (5),  and 

H-^^i=V,.(r)-V-=i=  (6) 

2"*VPi(r)  VPiW 

The  Euler-Lagrange  equation  for  the  two-body  correlations  can  be  formulated  as  an  equa¬ 
tion  for  the  inverse  of  the  static  structure  function  5(ri,r2)  (c.f.  Eq.  (B.2)  of  I): 


[S  ‘  * /fi  *  S  ■‘](ri,r2)  =  2Kp_fc(ri,r2)  +  [5;:‘  * *  Sf  *](ri,r2) 


is  the  particle-hole  interaction.  It  can  be  expressed  in  terms  of  the  Coulomb 
potential  Vcdri  -  r2{)  and  the  diagrammatic  quantities  defined  in  the  FHNC  theory.  Its 
precise  analytic  form  may  be  found  in  I.  For  our  present  purposes  it  is  sufficient  to  write 
^p-/i(ri,r2)  in  the  form 


^P-H(ri,r2)  =  \/pi(ri)t;c(iri  -  r2i)[l  -  G(ri.r2)i\/pi(r2). 


The  form  (8)  introduces  the  static  screening  factor®  1  -  G(ri,r2).  We  find  it  for  inhomo¬ 
geneous  systems  most  convenient  to  formulate  screening  in  configuration  space. 

The  total  energy  in  the  FHNC  approximation  used  in  I  has  the  form 


E  -  Tp  ^  Ee^  El  +  Ec  =  Ehf  +  Ec, 


where  Tf  is  the  kinetic  energy  of  the  model  system,  Ees  is  the  electrostatic  energy,  and  Ez 
is  the  exchange  energy,  all  calculated  in  the  single-particle  basis  {0t}-  Ec  is  the  correlation 
energy.  We  have  expressed  the  correlation  energy  as 


Er  =  ^Er,ot  ^Tjp  ^  Erpa  =  Erpa  t  A£^c 


l7*  M 


where 


(11) 


Erpa  =  -  J  *  (S'*  -  s;i)  .  if,  .  (5-‘  -  5;')!(r,r). 

The  analytic  forms  for  AEpot  and  ATjp  are  given  in  I;  they  are  irrelevant  for  our 
present  purpose.  Erpj^  is  the  RPA  contribution  to  the  correlation  energy  calculated 
with  the  screened  particle-hole  interaction  (8).  Setting  G{ri,r2)  =  0  in  Eq.  (8)  and 
ignoring  AEc  corresponds  to  a  “collective”  RPA,  which  is  equivalent  to  the  ordinary  RPA 
if  the  particle-hole  spectrum  is  replaced  by  an  effective  collective  mode^.  This  “collective” 
approximation  reproduces,  in  the  bulk  system,  the  RPA  for  the  correlation  energy  and  the 
static  structure  function  within  a  few  precent*.  In  the  RPA,  Eq.  (7)  determines  the  static 
structure  function  5(ri,r2),  which  is  then  used  to  calculate  the  correlation  energy  via  Eq. 
(11). 

In  I,  we  have  chosen  to  consider  jeliium  slabs  which  are  translationally  invariant  in 
the  X  -  y  plane  and  symmetric  about  2  =  0.  We  have  solved  the  coupled  FHNC/Euler- 
Lagrange  -  Hartree-Fock  equations  for  values  of  r,  corresponding  to  the  metals  At,  Pb,  Mg, 
Li.  .\a.  and  K  and  for  slab  dimensions  d  =  8,  10,  12.  and  14  oor,.  (an  is  the  Bohr  radius). 
From  the  ground  state  energy  as  a  function  of  the  particle  number  n  we  can  obtain  by  a 
linear  fit  20u  =  £<  (”)  -  Coo^  both  a  surface  energy  Ou  and  a  bulk  energy  per  particle,  Coo-  A 
corresponding  construction  can  be  made  independently  for  each  of  the  energy  contributions 
spelled  out  in  Eq.  (10). 

Table  1  shows,  for  the  above-mentioned  materials,  the  decomposition  of  the  total  jel¬ 
iium  surface  energy  into  contributions  originating  from  the  Hartree-Fock  energy  contri¬ 
bution  Ehf-  the  RP.A  energy  ErpaAorpa),  and  the  correlation  energy  correction 


A^c)(<7c)'  We  see  that  the  Hartree-Fock  ajid  the  RPA  contributions  dominate  Oc  by  an 
order  of  magnitude.  We  shall  therefore  ignore  <Jc  in  the  following  discussion. 

The  total  surface  energies  shown  in  Table  1  are  throughout  higher  than  those  obtained 
by  Lang  and  Kohn*.  In  order  to  understand  the  difference,  we  have  performed  an  array  of 
test-calculations  using  different  single-particle  wave  functions  and  energy  approximations 
(F«g-  1): 

(i)  In  order  to  test  the  dependence  of  the  FHNC  approximation  for  the  surface  energy  on 
the  single-particle  wave  functions,  we  have  performed  optimized  FHNC  calculations  with 
single-particle  wave  functions  generated  by  a  local  density-  functional  calculation  in  the 
sense  of  Lang  and  Kohn  from  the  bulk  FHNC  equation  of  state.  We  found  that  the  surface 
energies  from  these  calculations  agree,  within  a  few  percent,  with  the  ones  obtained  in  I. 

(ii)  We  have  obtained  a  density-dependent  static  screening  function  Gbuik{p,  iri  -r2|)  from 
an  FHNC  calculation  for  the  bulk  electron  gas,  and  used  it  in  local  density  approximation 

<''(ri.r2)  w  G’ft„<fc(p(rcm),!ri  -  r2l),  (12) 

to  solve  Eq.  (7)  and  to  calculate  the  RPA  energy  (11)  (LDS  approximation).  Tcm  = 
(fi  -r2)  2  is  the  center  of  m«i.ss  coordinate  of  the  two  particles.  To  guarantee  the  agreement 
between  the  kinetic,  Coulomb,  and  exchange  energies,  we  have  used  the  single-particle 
orbitals  obtained  in  (i).  We  found  that  the  surface  energies  obtained  in  this  way  are 
significantly  lower  than  the  FHNC  results  of  I,  but  in  good  agreement  with  the  Lang-Kohn 
results. 

(iii)  We  have  calculated,  in  our  geometry,  the  surface  energy  by  the  Lang-Kohn  procedure, 
using  the  bulk  FHNC  and  RP.-\  energy  per  particle.  We  found  surface  energies  in  agreement 


with  the  Lang-Kohn  results  within  a  few  percent. 

(iv)  We  have  performed  a  variational  RPA  calculation,  i.e.  solved  equation  (7)  with 
G'{ri,r2)  =  0  and  obtained  the  correlation  energy  from  Eq.  (11).  We  found  that  the 
surface  energies  agree  within  a  few  percent  with  the  LDS  calculation  (ii)  and  the  LDA 

(iii) . 

The  agreement  between  the  surface  energy  obtained  in  the  LDA  from  different  equa¬ 
tions  of  state,  as  described  in  (iii)  shows  that  the  surface  energy  does  not  depend  sensitively 
on  the  equation  of  state.  The  agreement  described  in  (i)  similarly  shows  that  the  correla¬ 
tion  contribution  to  the  surface  energy  is  insensitive  to  the  details  of  the  one-body  density. 

The  excellent  agreement  between  the  surface  energy  obtained  in  the  RPA  (procedure 

(iv) )  and  the  LDA  is  unexpected  and  at  present  unexplored.  Our  result  differs  somewhat 
from  the  conclusions  drawn  from  studies  of  surface-plasmon  contributions®’ to  the  surface 
energy.  The  agreement  deserves  further  investigation,  but  it  is  not  the  main  point  of  our 
study. 

('alculalion  (ii)  finally  traces  the  source  of  the  difference  between  the  full  FHNC 
calculation  of  1.  and  the  LDA:  This  difference  is  caused,  apart  from  the  small  correction 
<7,.  by  the  use  of  different  static  screening  functions  G(rj,r2).  We  conclude  that  this 
func'ion  (and  the  related  particle-hole  interaction  K,,_fc(ri.r2)  is  poorly  described  by  the 
local  density  approximation  (12).  The  difference  between  the  static  screening  function 
in  the  LDS  and  in  the  inhomogeneous  FHNC  is  not  very  large:  We  show  in  Fig.  2,  for 
\(j.  the  local  screening  function  obtained  from  the  inhomogeneous  FHNC,  parallel  and 
perpendicular  to  the  surface,  and  the  LDA  (12)  for  the  same  function,  for  different  values 
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of  the  center  of  mass  coordinate  2cm-  The  three  sets  of  curves  correspond  to  the  center  of 


mass  located  at  0.41aor3  outside  the  jeilium  edge,  on  the  jellium  edge,  and  O.Slaora  inside 


the  jellium  edge.  The  last  value  corresponds  to  the  maximum  of  the  first  Friedel  oscillation. 


The  basic  observation  is  that  G(ri,r2),  obtained  in  the  inhomogeneous  FHNC,  is  closer  to 


its  form  in  the  bulk  as  the  LDA.  The  rather  small  difference  in  the  local  screening  function 


induces,  however,  a  significant  change  of  the  surface  energy. 


We  come,  therefore,  to  the  main  conclusion  ot  this  work:  The  local  density  approx¬ 


imation  for  the  particle-hole  interaction  is  inadequate  in  metal  surfaces.  While  we  have 


confidence  in  the  qualitative  significance  of  our  observation  we  are,  due  to  the  delicacy 


of  the  effect,  presently  not  prepared  to  estimate  the  quantitative  accuracy  of  our  results. 


The  FHNC//0  approximation  used  here  is  the  minumum  level  at  which  the  effect  can  be 


studied,  but  it  is  improvable.  It  is  also  not  clear  whether  the  assumption  of  local  and  static 


screening  is  sufficient. 


Our  conclusion  that,  the  LDA  for  the  particle-hole  interaction  is  inadequate  can  be 


made  quite  plausible  sonsidering  the  surface  of  a  self-bound  system  like  liquid  ^He  or  ’He. 


.A  homogeneous  phase  of  these  systems  does  not  exist  below  the  density  where  the  velocity 


of  sound  vanishes.  .Any  LD.A  for  the  particle-hole  interaction  that  one  might  invent  must 


therefore  predict  an  instability  against  density  fluctuations.  Nevertheless,  the  particle- 


h<^le  interaction  of  two  particles  in  a  surface  should  never  predict  an  instability,  it  must 


therefore  deviate  significantly  from  any  LDA.  We  have  obviously  seen  that  this  effect,  while 


not  mandatory  for  stability,  is  also  significant  in  the  electron  liquid  surface. 


We  have  identified  the  point  where  the  LDS  fails  for  metal  surfaces,  and  what  needs  to 
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be  done  to  improve  upon  the  LDA  in  these  systems.  Simplifications  of  the  FHNC  theory 
may  now  be  found  a  posteriori  that  would  allow  the  same  calculational  accuracy  with 
simpler  computational  tools. 
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Metal 

OHF 

orpa 

Oc 

Ou 

2.07 

A1 

-1273 

994 

57 

-222 

2.30 

Pb 

-674 

812 

43 

181 

2.66 

Mg 

-215 

561 

37 

383 

3.28 

Li 

5 

328 

27 

360 

3.99 

Na 

40 

200 

21 

261 

4.96 

K 

35 

116 

8 

159 

TABLE  1 


Decomposition  of  the  total  surface  energy  into  its  contribution  from  the  Hartree- 
Fork  energy  ohf^  ’he  RPA  part,  of  the  total  energy  orpa-,  and  the  correction  <7c. 


Figure  Captions 


Fig.  1  The  surface  energy  of  jellium  slabs  obtained  in  the  different  procedures  described 
in  the  text  is  shown  as  a  function  of  r^;  The  FHNC  result  of  I  (solid  line), 
the  RPA  with  the  LDA  (12)  for  the  local  screening  function  (long-dashed  line), 
the  Lang-Kohn  result  (dash-dotted  line),  and  the  variational  RPA  (short-dashed 
line). 

Fig.  2  The  local  screening  function  G(ri,r2)  is  shown,  at  different  values  of  the  center  of 
mass-coordinate  2cm  relative  to  the  jellium  edge,  parallel  (solid  line)  and  perpen¬ 
dicular  (dashed  line)  to  the  surface.  Also  shown  is  the  local  density  approximation 
for  the  same  function  (circles). 
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